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Abstract 

We consider the L-state cyclic solid-on-solid lattice models under a class of 
open boundary conditions. The integrable boundary face weights are obtained 
by solving the reflection equations. Functional relations for the fused transfer 
matrices are presented for both periodic and open boundary conditions. The eigen- 
spectra of the unfused transfer matrix is obtained from the functional relations 
using the analytic Bethe ansatz. For a special case of crossing parameter A = vr/L, 
the finite-size corrections to the eigen-spectra of the critical models are obtained, 
from which the corresponding conformal dimensions follow. The calculation of the 
surface free energy away from criticality yields two surface specific heat exponents, 
Os = 2 — L/2£ and ai = 1 — L/l, where I = 1, 2, ■ ■ ■ , L — 1 coprime to L. These 
results are in agreement with the scaling relations as = on, + v and ai = a;, — 1. 
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1 Introduction 



Square lattice models in statistical mechanics with non-periodic boundary conditions have 
received intermittent attention over the years (see, e.g. [1-6]). Up until quite recently 
the systematic study of the integrability of such non-periodic systems lagged well behind 
the study of the corresponding periodic systems. It is well understood that models with 
periodic boundary conditions are integrable when their bulk/Boltzmann weights satisfy 
the Yang-Baxter equation 0. Since Sklyanin's work we now understand that lattice 
models with open boundary conditions are integrable if in addition the boundary weights 
satisfy the reflection equations |p. In particular, Sklyanin formulated the construction 
of commuting transfer matrices for the six-vertex model with open boundary conditions, 
from which the integrability is assured. Recent lectures on subsequent developments can 
be found in 0, [ill. 



Beyond their intrinsic mathematical interest, exactly solvable lattice models with open 
boundary conditions are attractive from the viewpoint of studying various surface critical 
phenomena [12-21]. Our motivation here is to study the surface critical behaviour of 
square lattice cyclic solid-on-solid (CSOS) models ||2^, These models are face models 
in which the adjacency condition between neighbouring heights is defined by the Dynkin 
diagram of the affine ^^li algebra. The CSOS model has been well-studied for periodic 
boundary conditions. The free energy, the local height probabilities and the correlation 
length have all been evaluated, along with their corresponding bulk critical exponents [22- 
25]. The complete operator content has been discussed in [^ and the fusion procedure 



has been carried out in ^7\, Some surface properties have been derived in 



The crossing or anisotropy parameter of the CSOS models is defined by A = in/L, 



where £ = 1, 2, ■ ■ ■ , L — 1 is coprime to L [23]. A special case of interest is L = 3 and 



£ = 2 which is related to Baxter's three colourings of the square lattice [^, 0, |2^ . 

The layout of this paper is the following. In section 2 the CSOS models with both 
periodic and open boundary conditions are described. We solve the reflection equations 
for the boundary face weights. The functional relations of the fused transfer matrices are 
also presented. In section |^ the eigen-spectra of the unfused transfer matrix is extracted 
from the functional relations following the analytic Bethe ansatz method. The finite-size 
corrections to the transfer matrix eigen-spectra at criticality are obtained for a special 
value of the crossing parameter. In section ^ the free energy of the open boundary 
models is shown to satisfy a unitarity relation. We solve the unitarity relation following 
the inversion relation method [0, ^ . From the singular part of the free energy we obtain 
two surface specific heat exponents in agreement with scaling predictions. We conclude 
with a brief discussion. 

2 CSOS models 



The CSOS lattice models ||22|, |23[ are a family of L-state face models |^ built on the affine 
Ai^_'^ Dynkin diagram. States at adjacent sites of the square lattice must be adjacent 
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on the Dynkin diagram. The cychc nature of the heights distinguishes the CSOS model 
from the corresponding RSOS model ||3^ built on the Al Dynkin diagram. 



2.1 Bulk face weights 

The allowed, or non-zero, face weights of the CSOS models are given by p2|, p3 
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where Wa = aX + wq. The height a = 1, 2, ■ ■ ■ , L and < ifo < vr is a free parameter. 
The crossing parameter A is given hj X = in/L, where i = 1,2,---,L — 1 is coprime to 
L and L > 2. The elliptic functions ^i{u), '&a{u) are standard theta functions of nome p 
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where < p < 1 with p = at criticality. 

These face weights satisfy the star-triangle equation 
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inversion/unitarity relations 



d g 
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and the crossing unitarity relations 
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where p{u) = ^i{X - u)^i{X + u)/^l{X). 
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2.2 Periodic boundaries 



There is a hierarchy of commuting families of transfer matrices constructed by the 
fusion procedure on the CSOS models under periodic boundary conditions. Let 
Wmxn(^^ ^ be the fused face weights with fusion level m and n in the vertical and 

horizontal directions, respectively |27, Then the fused transfer matrices T^'^'^\u) 
are defined with elements 



N 



(2.7) 



with oat+i = Oi and 6Ar+i = hi where is the number of faces in a row of the lattice. By 
construction the fused face weights satisfy the star-triangle equations, resulting in the 
commutation relations 

[T("'")(u),T("'")(w)] =0. (2.8) 

These fused transfer matrices satisfy groups of functional relations, which can be easily 
proved by fusion. Let us define 



along with the function 
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Then for each m = 1,2, 



the functional relations are 
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The unfused models of interest here are recovered by setting the fusion level to n 
m = 1. Disregarding finite-size corrections, the bulk free energy in this case satisfies 



nu)T{u + X) = fl 
This is the unitarity relation for periodic boundary conditions. 



(2.12) 



2.3 Boundary face weights 

Integrable models with open boundary conditions are defined by both the bulk and the 
boundary face weights. The latter are represented by three heights interacting round a 



triangular face |11, 34, 35, 36, 37]. For the CSOS models. 
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which satisfy the boundary version of the star-triangle equation (reflection equations) 



f,9 



a b 
9 c 



u — vjKyg 



f 

a b 
9 f 



u;nw 



a g 
d f 



u + v ]K[g 



u 



f 



+ v]K{d 



f 



u-nw 



a g 
d e 



u 



(2.14) 



In general there may be more arbitary parameters than ^. Inserting the CSOS bulk face 
weights ( |2.2| ) into the reflection equations and making use of the elliptic function identity 



(2.15) 



'i9i(x + y)di{x — y)'di{w + v)'di{w — v) — diiv + y)di{v — y)'di{w + x)'d4^{w — x) 
= 'di{x + v)'di{x — v)'d/i{w + y)'dii{w — y) 

we find the following CSOS boundary face weights 
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The identity ( p. 151) also plays a role in establishing the integrability of the bulk weights 
23[. 



It is obvious that the boundary face weights satisfy the crossing symmetry 
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2.4 Fusion results 



The fused transfer matrices T*^™''"^('u) of the open boundary CSOS models are defined by 
the following elements 
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where the right boundary fused face weights K^'' are given by fusing 
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while the left boundary fused face weig hts K^+^ are given by fusing 

t 
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-u + A; ^ 



These fused transfer matrices form commuting families, 







(2.22) 
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Like its periodic counterpart, the unfused transfer matrix is recovered by setting the 
fusion levels to m = n = 1, i.e., T{u) = T^^'^\u). 

The fusion procedure for face models with open boundary conditions has been demon- 
strated elsewhere ^ and thus we do not repeat the details here. However, it is 
worthwhile to write down the functional relations satisfied by the fused CSOS transfer 
matrices. We find 
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The boundaries contribute the diagonal matrix factors uj~{u) and uj^{u), with 
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Here height Cr (q) is located on the right (left) boundary. The matrix functions u;^(m) 
are simplified under the crossing symmetry ( p.l9| ), with 

^i{2X±2u) 
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Then the functional relations have similar forms to the periodic boundary case, namely 
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In this way, again after dropping the finite-size corrections, we arrive at the unitarity 
relation 

T{u)T{u + X) = f{u) , (2.31) 

where the function / is given by f{u) = fliu) p{2u)'d\{X) / 'd\{2X) after appropriate renor- 
malization of the free energies. This relation is sufficient to determine both the bulk and 
surface free energies |]37 |. 
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3 Exact Bethe ansatz solution 



Sklyanin presented the algebraic Bethe ansatz solution of the six- vertex model or spin- 
I XXZ chain with open boundary conditions. Unfortunately, the generalization of the 
algebraic Bethe ansatz to treat other integrable open boundary models has not made 
much progress, in particular, for models in which the arrow or spin reversal symmetry 
is broken. However, when such symmetry holds, the Bethe ansatz solutions of many 
integrable open boundary models have been obtained (see, e.g. [1^, Q). Here we 



show that the analytic ansatz method |^ can be applied to find the transfer matrix 
eigen-spectra of the CSOS models. This method has also been applied with success to 
the CSOS models with periodic boundary conditions pi|. 



3.1 Bethe ansatz solution 

Consider the transfer matrix T 



u 



which is the transfer matrix T^^'^^ with 



oo a]\f 
bo bN 

fixed heights bN,aN along the right boundary and bo,ao along the left boundary. For 
the solutions ( |2.16|) -( p.l8|) the transfer matrix is nonzero only for bo = ao and b^ = a^. 

Suppose Tj,^ii{u) = 
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for the eigen spectra of the transfer matrix Tb^d{u). We can check that the above ansatz 
satisfies the functional relation ( |2.3CI| ) if the parameters Uj satisfy 



TbAuk) = k = l,2, 



M. 



(3.2) 



As a result, the ansatz ( |3.1| ) should give the eigen-spectra of the CSOS transfer matrix, 
with (^.21) as the related Bethe ansatz equations. The eigen-spectra of the fused transfer 



matrices can been written down according to the fusion results from dOP and (3.2) 



3.2 Finite-size corrections 

At criticality the CSOS models are equivalent to the six-vertex model. The transfer 
matrix Tf,^d{u) is then independent of the boundary heights (as are the fused transfer 
matrices) and the Bethe ansatz solution (|3.1j)-(P^ is of the same form as the solution 
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of the open boundary six- vertex model [§, |T4|, |T^. Now the finite-size corrections to 



the fused transfer matrices of the f/q(s/2)-invariant six-vertex model have been derived 
T^|. Thus by a similar calculation we obtain the finite-size corrections to the fused 



m 



transfer matrices of the critical CSOS models. 

For simplicity we consider the limit C,± — ^ ^oo with A = vr/L and even. The fused 
CSOS transfer matrix eigenvalues then behave like 

logT(P'P)(«) = -2iV/,(«)-/,(«) + ^(c-24Ai,,,,)sin(Lw) + o(^). (3.3) 

The central charge is given by 

3p Qp 
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where p = 1,2, — 2 labels the fusion level. The conformal weights are given by 

_ {L - {L - p)rf - p^ z/(p-z/) 
^^'"'^ ALpiL-p) ~^2pip + 2) 

with 1/ a unique integer determined by 

r — 1 

u = r — 1 — p\ J . (3.6) 

p 

and r = 1, 3, ■ ■ ■ < L — 2. Here [xj is the largest integer part less than or equal to x. The 
functions fb{u) and fs{u) are, respectively, the bulk free energy and surface free energy 
of the critical models, which are also calculated for the off-critical models in the next 
section. They are not given explicitly here. 



4 Surface free energy and critical exponents 



The bulk and surface free energies of the CSOS models can be found from the unitarity 
relation ( p.31| ) with certain analyticity assumptions, as has been shown in the study of 
the eight- vertex model 0, |2g. 

The unitarity relation ( |2.31| ) combines the inversion relation and crossing symmetries 
of the local bulk and boundary face weights. We can separate the contributions from 
the bulk and surface free energies in this relation Let T{u) = Tb{u)Ts{u) be the 



^2N 



eigenvalues of the transfer matrix T{u). Define Tf, 
energies are defined by fb{u) = — log /«(,(«) and fs{u) = — log Kg 



and Tg = Ks, then the free 
u). We have 
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for the surface. Here height d [b) is located on the right (left) boundary. 

To solve the unitarity relations it is convenient to introduce the new variables 



X = e-^^l\ w 



-2™A^ q = e— 'A 



P 



along with the conjugate modulus transformation of the theta functions, 
The factor p{u, e) is harmless and will be disregarded, while 

oo 
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Suppose that Kb{w) is analytic and nonzero in the annulus < w < 1, we can Laurent 
expand fb{w) as logKf,(w) = I]J^L-oo Cn""^"- Then inserting the series expansion into the 



logarithm of both sides of (4J.) and equating coefficients in powers of w gives |23 
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Similarly, taking the Laurent expansion \og Hgiw) = CnW^ and solving the unitar- 

ity relation (Of) yields 
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The surface free energy is explicitly dependent on the boundary heights b, d and 

The specific heat critical exponents may be obtained from the leading order singularity 
of the free energies. In practice, the singular behaviour is extracted by means of the 



Poisson summation formula [0] . For the bulk free energy it follows that |]22|, p3 

fb{w,p) ^ P""^^ as p^O. 



(4.9) 



When i = 1 and L is even there is a multiplicative logp factor. It follows that the bulk 
specific heat exponent of the CSOS models is = 2 — ^ where we recall that \ = in/L. 
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The same idea can be applied to find the surface specific heat exponents from the surface 
energy fs{w,C,±,p)- Following |^0|, ^ and the treatment of the eight- vertex model 
we define the excess internal energy Cg, 



+ ei(p), 



(4.10) 



where ei {p) is called as the local internal energy in the surface layer, which is the correc- 
tion energy to the surface internal energy es{p), given by 



ei(p) 



dfs{w,^±,p) 



(4.11) 



The free parameter ^± appearing in the boundary face weights can be interpretted as a 
surface coupling. This has been explicitly shown in the study of the eight-vertex model 
pO| . The corresponding specific heats are defined by 



<9e,c 



dei 



dp ' dp 

Application of the Poisson summation formula to fs{w,$,±,p) yields 

es(p) ~ P^~^ 



(4.12) 

(4.13) 
(4.14) 



as p — s> 0, with a similar logp correction factor as for the bulk case. The surface specific 
heat exponents of the CSOS models follow as 
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For the three colouring problem (L = 3 and i = 2) we thus obtain the values = | and 



Recalling the bulk exponent at = 2 — ^ ^ and the correlation length exponent 



2A 



24] we are thus able to provide a further significant test of the scaling relations 



40, 41 



a. 



ai, + V and ai = a{, — 1 



(4.16) 



5 Discussion 

In this paper we have derived exact results for the critical surface properties of the CSOS 
lattice models. They can be generalized to the fused CSOS models. In this and related 
work HI the crossing unitarity relation plays a key role in deriving the surface free 

energy away from criticality. The CSOS lattice models with fixed boundary conditions 
on the square lattice with diagonal orientation can be treated in a similar manner by 
incorporating inhomogeneities into the bulk face weights and taking appropriate values of 
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the inhomogeneities and the boundary couphngs as has been exphcitly demonstrated 
for the ABF RSOS models However, we do not pursue this direction here as the 



change in lattice orientation does not effect the critical exponents. 

The excess surface critical exponent has been obtained from the singular leading 
term of the excess internal energy e^. It turns out that the singular leading term does 
not depend on the boundary face weights and thus as is independent of the details of 
the boundary weights. This behaviour has been already seen in the study of ABF model 



Other models of immediate interest are the dilute Al models [42, 43] which can be 
obtained from Kuniba's A2 face model |Q under restriction. The boundary face weights 



have been found for these models and the surface critical properties can thus be studied 



in a similar way EH . 
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